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Abstract 

Carrier interactions on graphene are studied. The study shows that besides the well known Coulomb 
repulsion between carriers, there also exist four-fermion interactions associated with U-process, one of which 
attracts carriers in different valleys. We then calculate the contributions to valley magnetic moment from 
vertex correction and from four-fermion corrections explicitly. The relative contributions are -18% and 3% 
respectively. At last we point out that we can mimic heavy quarkonium system by carrier interactions in 
graphene. 



1 Introduction 

Graphene [I], newly discovered two-dimensional 
crystals, has attracted more and more attentions 
of theorists and experimentalists "^ '^ . In graphene, 
there is a typical valley degeneracy, corresponding to 
the presence of two different valleys in the band struc- 
ture. However, as stated in the reference [B], such 
degeneracy makes it difficult to observe the intrinsic 
physics of a single valley in experiments How 
to distinguish carriers in the two valleys is therefore 
always a topic attracting literatures EEE] 

Ref. [9] pointed out that in close analogy with 
the spin degree, there is an intrinsic magnetic mo- 
ment associated with the valley index, which was 
called as valley magnetic moment (VMM). At tree 
level the valley magnetic moment is about 30 times 
that of the usual spin magnetic moment, therefore, 
"valleytronics" provides a new and much more stan- 
dard pathway to potential applications in a broad 
class of semiconductors as compared with the novel 
valley device in graphene nanoribbon However, 
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since in graphene the effective coupling ~ 1, a 
question to be posed is to what extent the calcula- 
tion in ref. ^ is valid. 

To answer the question, we first study carrier in- 
teractions. The study shows that, in tight binding 
approximation, besides the well-known Coulomb re- 
pulsion between electrons, there are also four-fermion 
interactions associated with U-process. The four- 
fermion interactions are type dependent and more 
significant, one of them attracts carriers in different 
valleys. Armed with the understanding of the in- 
teractions, we point out that there are two correc- 
tions to VMM at one-loop level. One is the vertex 
correction and the other is the four-fermion correc- 
tion. The vertex correction is similar to the anoma- 
lous magnetic moment of a particle in quantum elec- 



trodynamics (QED) except that carrier interactions 
on graphene are not "Lorentz covariant". Therefore, 
such correction always appears even in a one-valley 
system. Meanwhile, the valley degree is similar to 
the flavor degree in particle physics or high-energy 
physics. To compute anomalous magnetic moment 
of a particle due to flavor degree, one should also 
consider the weak interaction, an interaction between 
flavor degree. Our Yakawa-like four-fermion interac- 
tions are similar to the lower-energy approximation 
of the weak interactions. In this way, the correction 
to VMM due to valley degree appears at one-loop 
level. In contrast, such correction can not occur in 
QED. 

Our study shows that the total correction is about 
— 15%. Furthermore, since the corrections are inde- 
pendent on the divergence of the loop calculations, 
VMM can be used to check the validity of the per- 
turbational calculation. 



2 Carrier interactions 

Here we study carrier interactions. The study 
shows that besides the well known Coulomb repul- 
sion, there are also four-fermion interactions between 
carriers at different valleys, which are not only short- 
range but also contacting ones. 

For simplicity, we set A = 1 and X(r — r^/) the 
normalized orbital wave function of electron bound 
to atom A' , i.e. it satisfies / dvX{r — rA')X{v — 
ta") = &A'A" A-electron wave function ■i/'A(k) 



in position space is i/'A(k) = ^/u)/{2-k) 



,ikr- 



ta), where to = \/?ia? /2 is the area of the hexago- 
nal cell. For B-electron the case is similar. We then 
have < -(/)A;,(k')|V'Ao(k) >= (5AoA{,5(k - k'), where 



1 



To study carrier interactions, we consider 



where we have inserted the effective permittivity e in 
the last equation to include screening effect. We thus 
^(k) = ■^^7^'0A;(k2+k)V'A2(k2)t^'i/'A', (ki-k)?/>Ai(ki)pbtain the weU known Coulomb interaction. The 

type of A2 does not influence the results, that is. 



(2^)^ 



where A,, A{, A2, A'^ equal to A or B. If we ig- ^^6 Coulomb repulsion works both for carriers in the 

, ., ,. . .-v.. same valley and for carriers in the different valleys, 

nore mterchange mteractions, the mam contribution 

to X* (r' - v^, )X(r' - va,)VX' (r - v^, )X{v - va, ) "^^^^ Coulomb interactions, 

, , , , , ^. . ., ^ , there are other interactions which is related to valley- 

should be at vicinity r^' = ta, , r^' =rA, ,r « ta, , ... o , • . ■• 1 . 

and r ~ r We et ^ valley transition, buch interactions correspond to a 

^' U-process and therefore k ~ ±(47r/3a,0). To deal 

) Loe^ with such case, we substitute k-|- (Ij, 0) for k in Eq. 

iV|rAi -rA,|©- 

We first consider A-B interactions, that is, A2 = 



V{1^) = Sa,a', 



>5, 



- 5aiA'^5a^a'^ 



k{r 



,(2) 



B in Eq. ©. We get then 



where we have fixed A\ in the last step. We shall 
ignore the two delta functions thereinafter. 




|a2|. 
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K(k) 
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dr- 



27re^ 
ek 



(3a) 



, 2a ■sr^ 



^ifcy\/3a(n2 — ■!■) 



P2 



(3b) 



Compared to the long- wavelength result in Eq. (I3a|l , 
the valley-valley interaction suffers a coefficient sup- 
pression due to the large momentum transfer. How- 
ever, since such valley-valley interactions are short- 
range, it is not needed to consider screening effect. 
We therefore does not insert the effective permittiv- 
ity e in the above equation. 

Whereas when A2 = A, one should subtract 
the contribution from self interactions, which corre- 
sponds to (ni,n2) = (0,0) in Eq. ([2)). The result is 
then 



Figure 1: Graphene hexagonal lattice constructed 
as a superposition of two triangular lattices A and 
B, with bases vectors ai = a(i, -^) and a2 = 

a(^, — -^), where lattice constant a = |ai| = 
The reciprocal lattice vectors are bi ~ ^(1 

and b2 = ^(1, — -^) respectively. 

Without loss of generality, we set Ai — A. 
compute interactions between carriers, we first mark 
coordinates of A and B with two integers ni and 
712. From Fig[l] the coordinate of one atom A is 
set as (0,0). Then, for infinitely large graphene, co- 
ordinates of atom A are depicted as (ni/2, \/3(?ii — 
2n2)/2)a and coordinates of atoms B (ni/2, \/^{ni — 
2n2 + |)/2)a respectively, where a is lattice constant 
and ni, 712 are arbitrary integers. 

We put our focus on the interactions between elec- 
trons around ±K = ±(47r/3a,0). We first study the 
case where there is no valley-valley transition dur- 
ing interactions. For this case, we suppose |ka| ^ 1. 
In Eq. ([2]) the function in the summation is a slow- 
moving function, therefore, the summation can be 
replaced by an integral, 



K(k) 



-1.55e a. 



(3c) 



-1.55 is due to 



Here, the large negative coefficient 
the subtraction. 

Since Coulomb interaction is long-range, it does 
not depend on the distributing detail of the adjoint 
electrons. Thus, as shown in Eq. (|3a|l . such interac- 
tion is type-independent. However, the four-fermion 
form of the U-process implies that such interactions 
are short-range and they therefore depend on the dis- 
tributing detail. Therefore, as shown by (|3bp and 
(|3cp . such interactions are type-dependent. Refer- 
ence [13]!14! also proposed four-fermion interactions 
from different aspects. In Ref. [13| the authors add a 
near-neighbor interaction term and then, when they 
carry out momentum integral in the first Brillouin 
zone, they adhere the short-range interaction with 
the usual Coulomb interaction at |k| = f-^^, where 
k is the transfer momentum. In contrast, in our ap- 
proach, there is no artificial adhering and the interac- 
tions due to the valley transition are shown explicitly. 
Furthermore, the results obtained by our approach 
are suitable to take the quantum field theory calcu- 
lations. 

We emphasize that, besides the vertex correc- 
tion, Vs also contributes to VMM. Furthermore, since 
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Vs < 0, it takes attracting force between electrons in 
different valleys. The interaction rnay play crucial 
role in superconduction phenomena Therefore, 
Va deserves further research. 



3 The formal development of 
Lagrangian 

We first define two two-component spinors tp 

and V as follows: io — ( 'f^^^} ] and v = 

V ^k(p) J 



-i6_K(p) 



where ±K are two valleys. To de- 



ia_K(p) 

scribe the graphene dynamics in field theory lan- 
guage, we read the Lagrangian, 

£o = (pii'y^dt + iu7 • V — m):^ — e(^(7"j4" — jS'y ■ A) 
+x{W^dt -f- iu7 ■ V + m)x — ex(7"^" — Pi ■ ■ 

where /3 = w/c, w is the Fermi velocity of carriers, c is 

the effective light velocity in graphene, Ai — — {Vs — 

Vd)/2 and A2 = -{Vs + Vd)/2. We also set three 

/ 1 

gamma matrices as 7 = 70 = o"3 



into two pieces as follows: 

-Erlfril^A^^ — ■ Ar)iPr 

+ Xr{'^1^dt + iUr7 ■ V + mr)Xr 

-erXril'^A'^ - Prl ' Ar)Xr 

Air _ 1 _ A2r _ 2 - 

— — fr-y XrXrllH^r Y'^'''^ XrXrllVr 

-|-I^r(i'527°C^t + iSvVrJ ' V — 3m)(pr 

+Xr{iS2f°dt + iSyVr-y ■ V + 5m)Xr 

— er'fir{Sl'y°Ar — PrScf ■ Ar)ipr 

— erXr{Slj'Ar — /3r<5c7 ' A,.)Xr 

5lA _ 1 _ '52A _ 2 - ,r:\ 



= 5^71 = 



-71 = 7 (71 = 



. -i 
0' f 
-i 



and 
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7 =5 72 = 



72 = "/°a2 = y y ' ^^'^'^'^ 

ai 's are three Pauli matrices and metric matrix g^^ = 
diag{l, —1, —1} . Since four-fermion interactions in 
Eqs. H3b|) and pc|l are contacting ones, it is not nec- 
essary to introduce corresponding gauge field. In the 
above equation the energy gap m can be used to im- 
prove the use of graphene in making transistors and 
is therefore the one of the hot spots of literatures. In 
Ref. [16], the authors investigate the energy gap of 
graphene on a substrate BN, which is generated by 
the breaking of the A-B sublattice symmetry. How- 
ever, such energy gap has not been observed up to 
now. In Ref. [T7] the authors report that single layer 
graphene on SiC has a gap of 0.26eV, but the result 
is under debate [T8l[T9] . 

Utilizing the definitions of 7 matrix and 
boost generators corresponding to "Lorentz" trans- 
formation we find positive solution and nega- 
tive solution for Lp field (or negative solution 
and positive solution for x field) as u{p) = 
{^p° + m,v{p^ +ifi^)/\/p° +m)^ and v{p) = 
{v{p^ — ip^)/\/p'' -I- m, -^p" + m)"^ respectively. 
They meet ^u(p) = mu{p), v{p')0 = —v{p')m, where 
P = Ji^P'^ and p = up). 

However, the Lagrangian in Eq. Q is a bare one 
and it needs renormalization to match the observable 

quantities Having set (p — Z^'^pr^ X — 
1/2 

and A — Ar, where ipr, Xr and Ar are renormal- 
ized quantities, we split each term of the Lagrangian 



2fn — m,., 5^ — Z2v/vr — 1, Si — Z2Z^^^e/er — 1 = 
^- 1, = ZiP/Pr - 1, Six = ZiXi - Xir and 
— Z\\2 — \2r are determined by renormalized 
conditions. 

Since all the quantities in the following are renor- 
malized ones, all the subscripts r will be dropped 
out. 



4 Calculation of VMM 

To compute the VMM we first show Feynman 
rules in Fig. [2] (a)-(e). The contribution to VMM 
up to is depicted by Fig[2Ig), two diagrams 
in FiglUh), which is denoted by ^Ff, / = 1, 2, 
Fig[2ji), which is denoted by F^, and the countert- 
erm, Fig[2|f). The scattering amplitude of carrier 
under external gauge field is 



\M = -ie'{u(p')(l+5i)7^^(p)+r^+E'^rr}yl='(q), 

! = 1 

(6) 

where counterterm Si plays a similar role as Zint — 1 
and Z-kin — 1 in Ref. [21], and j4Jj'(q) is the Fourier 
transformation of the external field, p> and p are out- 
going and incoming momentums of carrier respec- 
tively. Since we are working in lower energy limit, 
we ignore the renormalization of Fermi velocity and 
charge. We get 



SVt = -iA, 



u{p')^\'$ - m)-y'^{ll:' - Tn)')iu{p) 
(27r)3 §2 



ir7][^'2 -m?] 



27rie^ 



ev^ J (27r)3 [^"2 _ + ir]][p - m2]|p - k| 

(7) 



respectively, where q — p' — p, k" — k + q and 
k' = k — q. In the above equation, we do not sum 
over the repeated index I and the terms proportional 
to ~ 10"* are neglected. 
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a b c 




Figure 2: Feynman rule and Feynman diagram 
on graphene. a) Propagator of Lp field with mo- 
mentum p, i (p)2^^'^_|_ig , where e is infinitesimal 
positive, b) Propagator of x field, i ^-^/j^+i^ - 
c) Propagator of gauge field, ^^zfj- d) Inter- 
action vertex between Lp field and gauge field, 
— ie'7^, where e' — e for /i = and e' — 
/3e for // = I, 2. For x fi^ld the interac- 
tion is similar, e) Two vertices of four-fermion 
interactions, and -i\2lf^i^l2j,j2- 

f)Counterterm vertex, — ie'5i7''. Since we are only 
concerned about the correction to VMM up to 
order e^, the renormalization of fermion velocity 
is ignored, g) Tree level diagram contributing to 
VMM. h) Four-fermion corrections, i) Vertex cor- 
rection. 

Since the external field is time-independent, q" = 

— p'' = in Eq. ([7]). If the electromagnetic field 
varies very slowly over a large region, Fourier compo- 
nents of the electromagnetic field will be concentrated 
about q = 0. We can thus take nonrelativistic limit, 
q — >■ 0, in iM, which means |pj, jp'|, |q| < m/v. 
Therefore, we have relations — (g)^ = v'^a^ > and 
W' ^ (py - v^p ■ p' ^ m\ 

To study the response to external magnetic field, 
we set time component of Ad as zero, i.e. Ac;(q) = 
(0, A(q)). We therefore only need to calculate the 
spatial part in iM. 

Since our theory violets the "Lorentz covari- 
ance", we should treat the result carefully. Fur- 
thermore, all the integrals in Eq. are divergent 
and therefore the result seems ambiguous. How- 
ever, we have the good news that the ambiguity 
have no effect on VMM. After a lengthy calcula- 
tion, such as Wick rotation ^ and the expansion of 
the result to order |p|, |p'|, |q|, we write F* and ^F; 
as u{C'iY + C2ie'"'q"'f3/2)u in nonrelativistic limit, 
where e^^ = -e^i = 1, = ^ g, Ci and C2 
depend on F^, 5r\ and SF^. For all the cases, Ci is 
divergent while C2 is finite. Together with the tree 
diagram and counterterm, {1+Si+Ci{r^)+Ci{5ri) + 



Ci{ST\))'u-y''u should be fixed to match renormal- 
ization conditions. Comparing with the Born ap- 
proximation for scattering from a potential of car- 
rier nearly p, q — s> 0, wc find that 1 -I- 5i -|- Ci(F') -|- 
C'l {Sr\)+Ci{Sr2) is just the electric charge of carrier, 
in units of e. Due to this, we set the renormalization 
condition as 1 + 5i + Ci{r) + Ci{Sr{) + CiiSri) = 1 
at p, q — >■ 0. This renormalization condition corre- 
sponds with the fact that the carrier at lower energy 
(p — 0) possesses unit (renormalized) charge e when 
scattered under external potential which varies very 
slowly. 

For finite term C2, we have 

Ignoring term proportional to p + p', which is the 
contribution of the operator p • A + A ■ p in the stan- 
dard kinetic energy term of nonrelativistic quantum 
mechanics, we rewrite uYu term as 

u{p )'y^u(p) — > u—u. (9) 

m 2 

We obtain, then, 

iM = -i2me^^e^(l - ^ + /^(A, + A2))B^ 
2 m -iev Anv^ 

(10) 

where = —\{q^ Aj.i — A].i) is magnetic field per- 
pendicular to graphene, ,^ — (1, 0)"^ is two-component 
spinor and C^^C = 1/2 = S3 indicates that electron 
pseudo-spin is 1/2. 

We interpret M as the Born approximation to the 
scattering of the electron from a potential. The po- 
tential is just that of a magnetic moment interaction, 
V(x) = -^^(K)B^ where 

3,_-, evB ,^ m ,^ , 

is the carrier VMM parallel to at K valley. For 
the hole, we get the same value with a necessary mi- 
nus sign. Similarly, for carrier at — K valley, VMM is 
also the same with a minus sign. Such phenomenon is 
known as the broken inversion symmetry in Ref. [9]. 

By recovering h, the leading term of VMM is 
ehvP/m, which is also obtained by ref. |9 . How- 
ever, besides the leading term, there are also other 
contributions to VMM. The relative contribution to 
VMM are 

mVs _ , 1 1.55s ma. , , 
~ 47rft2^2 -"(.-4+ 4^^„ )' (. ) 

where a — e^/{ehv) ~ 0.73 when e = 3. Substituting 

a = 2.46^, V ^ W'^cm/s into Eq. (O, we find the 
relative modifications to VMM due to vertex correc- 
tion and four-fermion interactions are about —18% 
and 3% respectively if we choose m — 0.26eV. 

It looks strange that it is not Vd but Vs which 
contributes to VMM. Such behavior stems from the 
definition of x field. From the definition of ip field and 
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X field Vd only relates to interaction between carri- 
ers in different valley with the same pseudo-spin so 
that it does not contribute to VMM. On the contrary, 
Vs relates to interaction between carriers in different 
valley with the different pseudo-spin. Therefore, only 
Vs contributes to VMM. 



teractions in graphene. Therefore, the study on the 
heavy-quarkonium system can also be carried out in 
graphene. 

This work is supported by the Cultivation Fund of 
the Key Scientific and Technical Innovation Project- 
Ministry of Education of China (No. 708082). 



5 Discussions 

In this paper we have discussed the carrier in- 
teractions. The study reveals that besides the well 
known Coulomb repulsion between carriers, there are 
four-fermion interactions between carriers in different 
valleys. Since the interactions are short-range and 
contacting ones, they depend on the atom colloca- 
tion detail. Therefore, the four-fermion interactions 
are type dependent. Our study shows that one of the 
four-fermion interactions attracts carriers in different 
valleys, which we believe to be helpful in understand- 
ing the unusual superconduction effect in graphene. 

We also compute VMM from the tree level dia- 
gram, the vertex correction and the four-fermion in- 
teractions respectively. The contribution from the 
tree level diagram agrees with the result obtained in 
ref. [9]. The other two contributions counteract each 
other and therefore the total contribution to VMM 
is about —15% if we choose m — 0.26el/ and e = 3. 

The very high accurate measurement of spin mag- 
netic moment is very important, both from theoreti- 
cal viewpoint and from practical one. Similarly, our 
result on VMM is also significant to valleytronics in 
graphene, especially to the future apparatus design 
based on valleytronics. Our study also points out 
that, in close analogy to the Zeeman split, the con- 
tribution to VMM induced by Vs is inherent, since it 
is independent on the energy gap m. In other words, 
to measure the magnetic moment induced by Vs, we 
may choose the substrate freely, although different 
substrate may generate different energy gap and dif- 
ferent effective permittivity. 

From Eq. (|12p . a plays the same role as the 

2 

fine structure constant in QED, = ~ 1/137. 
However, because a is about 100 times larger than 
Qe, it is hard to state that we mimic QED by car- 
rier interaction in graphene. Meanwhile, when we 
deal with problems dominated by quantum chro- 
modynamics(QCD), especially in heavy quarkonium, 
such as cc system and bb system, we always take 

2 

as = 1^, where Qs is the QCD coupling, as the es- 
timate of the effectiveness of perturbational expan- 
sion. (In many cases when we deal with such prob- 
lem we take an approach very similar to QED, up 
to an unimportant color factor.) At energy scale 
740Mev, as{7AQMev) « 0.73 Ri a I^l]. Noticing that 
the energy scale is close to the soft scale of cc and 
bb systems the dynamics of which is depicted 
by nonrelativistic QCD, we conclude that we can 
mimic the heavy quarkonium system by carrier in- 
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